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Abstract 

This paper examines the impact of two important geometrical parameters, namely the thickness and source/drain 
extensions on the performance of low doped p-type double lateral gate junctionless transistors (DGJLTs). The three 
dimensional Technology Computer-Aided Design simulation is implemented to calculate the characteristics of the devices 
with different thickness and source/drain extension and based on that, the parameters such as threshold voltage, 
transconductance and resistance in saturation region are analyzed. In addition, simulation results provide a physical 
explanation for the variation of device characteristics given by the variation of geometric parameters, mainly based on 
investigation of the electric field components and the carries density variation. It is shown that, the variation of the carrier 
density is the main factor which affects the characteristics of the device when the device's thickness is varied. However, the 
electric field is mainly responsible for variation of the characteristics when the source/drain extension is changed. 
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Introduction 

As the conventional planar metal-oxide-semiconductor field- 
effect transistors (MOSFETs) dimensions scale down to tens of 
nanometers, the formation of sharp source/ channel and channel/ 
drain junctions becomes crucial for suppressing the issues such as 
leakage current and short-channel effect (SCE). In nano scale 
devices, on the other hand, the formation of the ultra sharp 
junctions imposes severe challenges on doping techniques due to 
the difficulty to control the distribution of dopants at the 
metallurgical junctions and the intrinsic discreteness of the dopant 
itself Recently, a novel device called junctionless transistor (JET) 
has gain significant attention due to the advantages associated to 
the specific design of the device and simplification of the 
fabrication process [1,2,3]. From the structural point of view, 
the JETs are heavily doped gated resistors with narrow silicon 
body in multiple gate architecture. Operationally, the JETs are 
volume depleted in the off state (at zero gate bias) due to the 
effective work-function difference between the gate and the 
channel [4]. Conduction mechanism in JETs only occurs in the 
bulk of the channel. Besides, no surface accumulation layer is 
formed below the threshold voltage [5] . The freedom of JETs from 
different types of doping and associated gradients is considered as 
the main advantage of this design which simplifies the fabrication 
process. In addition, in subthreshold region device can obtain 
enhanced immunity from short channel effects in comparison to 
conventional inversion mode devices [6] . Performance comparison 
between bulk and SOI JETs [7,8], temperature dependence [9], 
scalability [10,11], ballistic transport [12], and analog and digital 



applications of JETs [13,14] have been investigated in the 
literature. 

Recently, the researchers of the present study have reported the 
fabrication and experimental characterization of low doped p-type 
junctionless transistors with single and double lateral gate(s) 
through the unconventional method of atomic force microscopy 
(AFM) nanolithography [15,16,17]. Simulations of carrier's 
transport and pinch off mechanism of the device have been 
presented in Refs. [18,19]. Computational studies are absolutely 
critical to predict the ultimate performance of the nanoscale 
devices and to impel forward future experiments or technological 
developments. As a result, in this work, the impact of the two 
important geometric parameters to the operations of the p-type 
double lateral gate junctionless transistors (DGJET) structure is 
comprehensively explored. First, the device simulation structures 
and models implemented in the present research are presented. 
Then, the operation and performance of the device in on and off 
state is briefly explained. Finally, the impact of geometric 
parameters of the device and particularly active layer thickness 
and source/ drain extensions on characteristics of the device, are 
investigated based on the variation of hole density and electric field 
components. 

Methods 

The structure of a typical DGJET used for the simulations is 
shown in Fig. 1. Parameters used for the device simulations are 
listed in Table 1. The simulated structures have a uniform low 
doping concentration (Boron 10^^) throughout all active regions of 
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Figure 1. Isometric view of simulated DGJLT. 

doi:10.1371/journal.pone.0095182.g001 

the device. The complete structure sits on a 145 nm ideal oxide. 
Different zones (Xj, Xn, and Xjn) correspond to the source 
extension, the channel and the drain extension, respectively are 
indicated for enhancement of the discussion and will be referred 
throughout the discussion. 

The models used for the simulations are calibrated against a 
long-channel experimentally demonstrated DGJLT, as described 
in [18,19]. We use a hydrodynamic model in Sentaurus software 
D-20 10.03 [20] as the platform for the 3-D TCAD simulation 
presented in this study. Besides the fundamental equations, the 
doping-dependent Masetti mobility model and the default 
S ho ckley— Read-Hall (SRH) recombination-generation model is 
activated. SRH recombination-generation is applied in order to 
consider the leakage current and recombination through deep 
defect levels in the gap [21]. 

Results and Discussion 

On and off state behavior of the device 

Fig. 2 shows the simulated transfer characteristics for the 
nominal DGJLT with 100 nm thickness and width for Vi3= — 
0.05 V and —1.0 V. Zones Xj, Xn, and Xni have the doping 
density of 10^^ cm ^, and the dimensions are considered as 2 |im, 
200 nm, and 2 |im, respectively. The nominal simulated device is 
considered analogous to the fabricated device in order to give the 
basic view of the device operation. It can be seen that the device is 
in on state at zero gate voltage with /^^Z/o^ ratio of 10^. The device 
is driven into off state by increasing the positive gate voltage. 
However, negative gate voltage is not able to increase the current 



Table 1. Simulated devices parameters. 





Parameters 


Value 


Device layer thickness 


20, 40, 100 nm 


Zones X| and Xm lengtli 


500 nm, 1 |im and 2 |im 


Zone X|| length 


200 nm 


Contact work function 


5.12 eV 


Gate voltage 


-2 V to +2 V 


Drain voltage 


-0.05 V to -1 V 
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significantly. Since the device is in on state and already near the 
flatband condition at zero gate voltage, it can be claimed that the 
device is principally working similar to the pinch-off transistors 
[22,23]. In the on state, the majority carriers moves in the volume 
of the channel and the depletion mechanism of the device due to 
the positive gate voltage applied to the lateral gates starts from the 
bottom corners of the channel at the Si/BOX interface, face to the 
lateral gates, and expands toward the center and top surface of the 
channel [18]. The electrostatic squeezing of the channel creates a 
large barrier along the holes transmission path. Increasing the 
drain voltage creates a depletion layer due to the sweeping of the 
carriers and at very high drain voltage, this depletion area along 
the electric field works as a barrier and prevents the electrical field 
from the propagating into the channel. This creates the current 
saturation, even for Vq = 0 V. 

Influence of thickness 

In the p-type DGJLTs, the conduction mechanism is through 
the current of majority carriers (holes) flowing in the bulk of the 
channel when the device is in the on state, and this volume as a 
conduction path would be gradually shrunk from bottom surface 
at Si/BOX interface to the top of the channel by increasing the 
positive gate voltage to turn the device off. This implies that the 
thickness has a strong influence on the characteristics and switch 
behavior of the DGJLTs. In order to investigation the effect of 
thickness on the behavior of the DGJLTs, the structures with 
thickness of 20 and 40 nm were simulated and the results 
compared to the nominal simulated device with 100 nm thickness. 
To avoid any quantum mechanical confinement effects, the 
thickness lower than 20 nm has not been investigated in this study. 
The dependence of transfer characteristics of DGJLT on thickness 
variation, when the drain is biased by Vd = — 1 .0 V is shown in 
Fig. 3. The inset shows the variation of on and off current with the 
thickness. As it is expected, the device with 20 nm thickness 
provides the highest /^^/Zo^ratio (—5x10^) and a small decrease of 
the on state drain current is obtained for thinner devices. 

The threshold voltage variation over the change of the thickness 
(dVth/dTsi~5 mV/nm) was negligible, mainly due to the 
suppression of random impurity fluctuation in these low doped 
devices [24]. The low doped channel in thin SOI devices has been 
implemented experimentally to minimize threshold voltage 
variations caused by random impurity fluctuation effects [25]. 
Another source of threshold voltage fluctuation is the scattering of 
source/drain doping atoms into the channel. This source of 




vg(v) 

Figure 2. Transfer characteristics simulations of DGJLT with 
100 nm width and thickness, nanowire length of 4.2 jam, and 
channel length of 200 nm. 

doi:10.1371/joumal.pone.0095182.g002 
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Figure 3. Transfer characteristics simulations of DGJLTs with 
three different thiclcness (20, 40, 100 nm). The inset plots the on 
(Vg = 0 V) and off (Vg = 2.0 V) current variation with thickness, at Vd= - 
1.0 V. 

doi:10.1371/journal.pone.0095182.g003 

threshold voltage fluctuation is eliminated in DGJLTs because of 
the absence of junctions and doping gradients [26] . Based on the 
characteristics of the devices, we try to extract some of the 
important parameters in order to describe the variation of the 
device operation with respect to the thickness. 

The transconductance (g^) of the DGJLTs with three different 
thicknesses as a function of gate voltage is shown in Fig. 4a. The 
inset shows the variation of the g^ as a function of thickness at 
three different gate voltages. The overall trend demonstrates the 
variation of all devices from on to off state. Sharp decrease of 
transconductance shows the effect of positive gate voltage to 
deplete the channel according to the pinch off mechanism. 
Moreover, the lower transconductance in thinner device is 
attributed to the scarceness of the carriers flowing through the 
body of channel. Accordingly, the on state current is decreased in 
the thinner device. 

A useful function for the extraction of the parameters in the 
conventional devices is the combination of the drain current and 
transconductance. This parameter (iD/gm^'^) known as Y(Vg) is 
conventionally used to extract parameters such as, the threshold 
voltage and low field mobility [27]. The lr>/gm^'^ variation versus 
gate voltage for devices with different thickness is plotted in Fig. 4b. 
Before investigating this function, it worth noting that, the 
operational principles of p-type DGJLTs, is different with that of 
gated resistors JLTs. In the gated resistor JLTs, below threshold 
the device is fully depleted. The conversion of the channel from 
full depletion which no conduction occurs to partial depletion 
which the bulk conduction starts can be achieved by increasing the 
gate voltage. Further increase in the gate voltage to the flatband 
voltage (Vfb) makes channel to be completely neutral and above 
flatband an accumulation layer can be created at the interface with 
the gate [1,5]. On the other hand, the DGJLTs are normally on 
device which forced to the off state by increasing the gate voltage 
and depleting the channel from the majority carriers. The 
variation of majority carriers density in the cross section of 
channel by increasing the gate voltage is shown in [18]. It is shown 
that, by small increasing of the gate voltage, the gates can just 
influence on the corner of the channel and it cannot affect the 
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Figure 4. Variation of transconductance (gm) (a), plot of the 
ratio of the drain current to the square root of the transcon- 
ductance (io/gm^'^) (b) and gm/lo versus gate voltage (c) for 
three different thicknesses (20, 40, and 100 nm) in saturation 
(Vd=-1.0 V). 
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current significantly, mainly due to the lateral gates design and 
lack of the gate oxide. As a result, in l-o/gm^'^ ratio, the numerator 
cannot be changed significantly; however the decrement of g^ in 
denominator is very sharp which is appeared as an increase of 1^/ 
gm^'^ in Fig. 4b. On the other hand, by increasing the gate voltage, 
the lateral gate influence on the current is more significant, while 
the transconductance variation is negligible. This behavior is 
shown with a graduate decrease in lD/gm°'^ for higher gate 
voltage. It should be noticed that, the lower value of I^/ gm^'^ in 
thinner devices might be due the lower carrier concentration and 
lower value of transconductance; however in thinner devices the 
current increment is dominant. 

A plot of gm/Io as a function of gate voltage when the devices 
with different thickness are biased in the saturation region (Vd = ~ 
1.0 V), is shown in Fig. 4c. This ratio indeed is a universal 
characteristic of all the transistors belonging to a same process and 
known as a measure of the transconductance generation efficiency, 
since it is a measure of the efficiency of the device to transform 
current into transconductance [28]. When compared to the 
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Figure 5. Hole density as a function of position along a horizontal cut line (a) and hole density along a vertical cut at y = 0 (b), for 
devices with different thicknesses of 20, 40, 100 nm. Vd = -1.0 V, Vg = +2 V. 
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Figure 6. Electric field as a function of position along a horizontal cut line (a) and Electric field along a vertical cut at y = 0 (b), for 
devices with different thicknesses of 20, 40, 100 nm. Vd = -1.0 V, Vg = +2 V. 
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Figure 7. Transfer characteristics (Id-Vg) (a) and drain conduc- 
tance as a function of drain voltage (b), for different zones X|/ 
X||| length at Vg = 0.5 V. Zones X, and Xm have the length of 500, 
1000, and 2000 nm. The length of zone X,, is constant at 200 nm. 
doi:10.1371/journal.pone.0095182.g007 

degradation of drive current and transconductance in thinner 
devices, the degradation in gm/Io is seen to be smaller and gm/lD 
is nearly identical for the different thickness. This trend is 
attributed to the fact that all devices have similar efiiciency to 
transform current into transconductance. In addition, gm/Io in 
JLTs is principally controlled by the body factor of the device 
which is equal to unity in these devices [14]. In the next section, in 
order to give a deeper physical perception about the effect of 
thickness variation on the device performance, we investigate the 
parameters which directly influence the output characteristics. 

Origins of characteristics change by thickness variation 

In the DGLJTs, the gradient of carrier concentration along the 
source/drain and channel is zero, hence no diffusion can take 
place and the current is dominated by the drift current. For 
semiconductors with holes as majority carriers, the drift current 
under an applied field is given by [21]: 



J = q\ipFE 



hence the current can be written as 



P(A)EdA = qfiP(x)E(x) 



(1) 



(2) 



where q is the electron charge, |i is mobility of carriers, P(x) is the 
carriers density (p(A)) integrated over the cross sectional area (A) of 




1.0 1.5 

Xj/Xjii (|Lim) 

Figure 8. Total resistance versus zone Xi/Xm length for two 
different gate voltages of 0.0 and 0.5 V, in the saturation 
region. 

doi:1 0.1 371/journal.pone.00951 82.g008 

the device, and E is the electric field. It can be seen that, mobility, 
carrier density and electric field are factors which influence on 
output characteristics of the device, e.g. Id and gm- 

In the DGJLT the conduction path locates near the centre of 
the channel instead of being confined into the surface channel like 
conventional MOSFETs. This allows for the holes to move 
through the silicon with bulk mobility. In addition, due to the 
fabrication method based on the AFM nanolithography the 
surface and the body of the upper silicon layer in SOI remains 
intact. This makes the possibility to obtain more bulk properties, 
such as mobility value close to the bulk mobility even at very thin 
devices [15]. Accordingly, we can claim that the mobility variation 
is not the main factor for variation of characteristics of the device 
when thickness is varied. In this matter, as the most influential 
factors, the carrier density and electric field are taken into account 
for further investigation about thickness variation. 

The holes density (majority carriers) variation, along a 
longitudinal cutline and in a cross section of devices with three 
different thickness (20, 40, 100 nm) at Vq = +2 V and Vd = - 1 V 
is presented in Fig. 5 (a-b), respectively. The effect of gate voltage 
variation on hole density for the nominal device was already 
investigated in Ref [18]. In this gate voltage, the lateral gates have 
the highest impact on the channel with maximum variation of the 
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Figure 9. Electric field along the horizontal cut line for different 
zones X| and Xm lengths, Vd= -1.0 V, Vg=+2 V. 

doi:10.1371/journal.pone.0095182.g009 
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Figure 10. Simulated hole density as a function of position along a cut at Z = 50 nm (a) and hole density as a function of position 
along a horizontal cut line (b), for different zones X| and Xm lengths at Vd= — 1.0 V and Vg = +2 V. 

doi:1 0.1 371 /journal.pone.00951 82.g01 0 



carrier population. It is observed that, as the channel thickness 
decreased, the lateral gates were able to vary the carriers more 
effectively and consequently the current variation is more 
significant in the thinner devices. This trend confirms the variation 
of gm results already presented in Fig. 4a, where lower g^ for 
smaller thickness can be explained by lower hole density in the 
channel. 

Figs. 6 (a-b) present the electric field along a horizontal line 
along the channel and in a cross section for devices with different 
thicknesses. It can be observed that, the peak of the electric field is 
occurred in the drain extension and not into the channel which 
confirms the typical behavior of the JLTs for all thicknesses [29] . 
The maximum electric field along the channel appeared in the 
device with the smallest thickness and for all devices the electric 
field reach to the lowest value at the centre of the channel (x = 0). 

In fact, the overall electric field behavior along the nanowire 
was consistent with the hole density. Wherever the electric field 
was stronger the higher rate of hole depletion occurred. This trend 
for device with the lower thickness was even more intensified. The 
stronger electric field and the highest hole depletion were achieved 
at the lowest thickness (20 nm). 

According to the definition of the transconductance and 

considering that the value of the mobility is almost constant, from 
the differentiation of Equation. 2, we can obtain: 

This shows that two factors can directly cause the transconduc- 
tance variation. The first factor is the variation of hole density with 
gate voltage multiplied by electric field along the channel 
(^y^ E(x)) and the second factor is the variation of electric field 
with gate voltage multiplied by hole density along the channel 
ffv^ P(x)). As it is presented in Figs 6a and 7 a the hole density is 

significantly reduced with decreasing the thickness, however the 
magnitude of electric field is increased with decreasing the 
thickness. As the value of g^ is lowered with decreasing the 
thickness (Fig. 4a), then it can be concluded that (^^^E{pc)) is the 

dominant factor contributing to the reduction of g^ in thinner 
devices. 

Influence of zones Xi/Xm length, Fixed Zone Xn 

According to the specific design of the device, the variation of 
source/drain extension lengths (Xj/Xni) is another important 
geometric parameter which can significantly influence on the 
output characteristics. Since the nominal device was experimen- 
tally characterized with the channel length (Xu) of 200 nm, the 
same Xn length is taken for all devices to purely investigate the 
effect of zones Xi/Xm length variation. The transfer character- 
istics (Id~Vg) and drain conductance variation as a function of 
drain voltage (gd-V^) of the devices with different zones Xj/Xjn 
length, in comparison with the nominal device, are presented in 
Figs. 7 (a-b). The device display typical characteristics over the 
investigated range of zones Xj and Xni lengths. On and o^currents 
were increased by decreasing the zones Xj and Xm lengths. 
According to the model proposed in [30], in addition to the 
channel cross section, the output current in the flatband condition 



depends on the effective drain voltage at the channel (Vch) and 
doping concentration. 

The transfer characteristics presented in Fig. 7aindicates that as 
the zones Xj/Xni length scaled down to 500 nm, the subthreshold 
leakage increased by about 2 orders of magnitude resulting in 
degradation of I^JI^ff ratio. However, the comparison of drain 
conductance as a function of drain voltage for different zones Xj/ 
Xxn length presented in Fig. 7b indicates that by decreasing zones 
Xj/Xni length an increasing trend of the drain conductance is 
observed for all drain voltages. In order to understand the reason 
for increase of the on and off state current when the zones Xj/Xni 
length decreases, some of the factors which influence on the 
characteristics of the device are analyzed. 

Total resistance (Rtot) versus zones Xj/Xni length obtained 
using the slope of the output characteristics data for devices with 
zones Xj/Xni length of 0.5 |im, 1 |im, and 2 |im at two different 
gate voltages are shown in Fig. 8. The total resistance is varied 
from 6 MQ for device with 0.5 |im Xi/Xm length to 28 MQ for 
device with 2 |im Xj/Xni length, at Vq = 0 V. The resistance 
obtained includes the nanowire resistance (Rnw)? contacts 
resistance (RJ and the resistance of the external circuit (Rext) [31], 

= + ^ext + ^Cl + ^C2 (4) 

The Rci and Rc2 could be estimated by simply extrapolating the 
curves to Z = 0 since Rnw decreases to zero and i^tot asymptot- 
ically approaches the value of R^ [32]. By decreasing the zone Xj/ 
Xjii the only part which varies is the Rnw due to the variation of 
length. Therefore, the lower Rnw the device with 0.5 |im result 
in the enhanced on -state performance. However, with decreasing 
length the off state current (I^^ is also increased, which is not 
favorable. The reason for I^^ value increment can be understood 
by investigating the electric field and hole density variation along 
the channel. 

The electric field as a function of position along the current 
direction in devices with three different zones Xi/Xm lengths is 
presented in Fig. 9. It is shown that, the interaction of electric fields 
from lateral gates and drain contact was diminished in the channel 
and enhanced in drain extension for all devices. However, for the 
device with the shortest length (Xj = Xni = 500 nm), the electric 
field was stronger at the drain extension compared to the devices 
with longer Xj/Xni length. Moreover, the peak of the electric field 
is located closer to the channel. For the short Xi/Xm length, the 
strong electric field from the drain could overcome the zone Xni 
and suck the holes directly from the channel or even zone Xj 



Table 2. Dependency of the devices' characteristics to the 
scaling of thickness and zones Xi/Xm length. 





Thickness scaling 


Zone X|/X||| scaling 


Threshold voltage (Vth) -Insensitive 


—Insensitive 


Leakage current Decreases 


Increases 


Hole density in the channel Decreases 


Increases 


Electric field Increases 


Increases 


The most Influential factor Hole density 


Electric field 



doi:10.1371/journal.pone.0095182.t002 
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(punching effect). The carriers which are dislocated by electric 
field need to be relocated somewhere in the channel (zone Xn) or 
in the zone Xni- This can vary the carrier configuration in 
different parts of the device. 

The holes density variation along the cut at Z = 50 nm and 
along a cutline as a function of position for the devices with three 
different zones Xj/Xni length are shown in Fig. 10 (a— b), 
respectively. It can be observed that the lowest depletion in the 
channel occurred for the shortest length (Xj = Xni = 500 nm). For 
the devices with longer length, the depleted carriers from channel 
could be spread into the zone Xjn- 

However, for the device with the shortest length, there is no 
area for the depleted holes to be located because of the punching 
effect and drain influence. This could provide the lower rate of 
depletion in the channel at off state and explain the higher 
subthreshold current for the device with the shortest length. 

Conclusion 

Thickness and source/ drain extension lengths as two important 
geometric parameters of the double lateral gate junctionless 
transistors (DGJLT) have been investigated using 3-D TCAD 
simulation tools. The device has shown a good scaling ability in the 
investigated parameters. Unlike the high doped JLTs, the 
threshold voltage remains approximately insensitive to the 
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